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Glosario 


amplitud vectorial del campo electromagnético; 
ángulo cenital (en radianes); 

ángulo azimutal (en radianes); 

fase; 

frecuencia temporal; 


función de onda vectorial del campo 
electromagnético; 


variación de la intensidad debida a la 
interferencia de N haces de luz; 

intensidad de la interferencia de dos haces m, n; 
longitud de onda; 

función de x, 


transformada de Fourier de la función x respecto 
ax; 

variable temporal (tiempo); 

variable temporal recíproca; 

variables espaciales cartesianas; 


variable espacial recíproca de x; 
variable espacial recíproca de y ; 


AS 


pa] 


IO 
IS 
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O, 


R”: 
AlfÍx]]: 
WLf lx]: 


r. 


P: 


vector del número de onda; 

vector en dirección radial; 

vector en dirección azimutal; 

vector en dirección cenital; 

vector unitario en dirección radial; 
vector unitario en dirección del ángulo 
azimutal; 

vector unitario en dirección del enésimo 
ángulo cenital; 

Espacio tridimensional (x, y, z); 
función de Ambigiedad; 


función de distribución de Wigner; 


variable espacial igual a ,/x” + ye 
variable espacial igual a A 


plano 7 —Ó : plano X-Y, en coordenadas polares 
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Introducción 


contexto actual y antecedentes del problema 
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La óptica es una de las ramas de la física mejor 
establecidas en términos del conjunto de 
conocimientos básicos requeridos para entender los 
fenómenos ópticos que aparecen en la Naturaleza. Es 
bien sabido que los primeros estudios de óptica datan 
de los antiguos griegos y que casi todos los grandes 
físicos de la historia trabajaron, ya sea directa o 
indirectamente, problemas ópticos. En este sentido, 
parecería que hay pocos principios novedosos qué 
estudiar en óptica, desde el punto de vista de la teoría 
fundamental. Sin embargo, los últimos años han visto 
renacer el interés de muchos grandes físicos 
contemporáneos por la óptica teórica, no solamente 
con métodos matemáticos más poderosos, sino con una 
serie de propuestas tecnológicas novedosas. 
Consideremos, por ejemplo, el bien conocido fenómeno 
de la interferencia, que es una de las bases de la 
llamada “óptica física". A pesar de que los fenómenos 
de interferencia de la luz se conocen desde Huygens, 
en el siglo XVI, y de que, en el siglo XIX, Young 
estudió con sumo detalle la interferencia constructiva 
y destructiva de haces luminosos, el principio de 
interferencia holográfica es relativamente reciente, 
ya que fue propuesto teóricamente por Gabor en 1948 
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[1] y los subsecuentes desarrollos [2] han abierto 
grandes perspectivas prácticas para la interferometría 
holográfica. En este contexto, el uso de más de dos 
haces que interfieran entre sí fue propuesta 
originalmente por Cowan en 1987 [3]-[7], obteniendo 
resultados experimentales con aplicaciones 
tecnológicas muy importantes, y no totalmente 
explotadas hasta la fecha. Análisis teóricos sobre 
interferencia múltiple se han reportado [8] y son bien 
conocidos, aunque el caso de simetría cilíndrica no ha 
sido presentado en la literatura especializada con 
detalle. En el caso de la difracción, fenómeno 
intimamente asociado a la interferencia, se han 
propuesto patrones de difracción basados en fractales 
[9] y análisis de Fourier, con la finalidad de hacer más 
eficiente la generación de patrones de interferencia 
[10], tanto a nivel teórico como tecnológico. La 
investigación tecnológica [11]-[26] ha abierto un 
mercado importante en la industria de la holografía 
tradicional y de la holografía digital y esto, entre 
otras razones, ha reavivado intensamente el actual 
interés en interferencia. 


Recientemente, se ha reportado el uso de técnicas de 
espacio fase, como la Transformada de Wigner, 


mediante funciones intermedias (como la 
Transformada de Fourier Fraccionaria) a problemas de 
Óptica [27]-[41] siendo de especial interés para esta 
tesis su uso en problemas de interferencia. La 
caracterización de sistemas ópticos como la 
interferencia múltiple, o la interferencia de anillos 
circulares puede comprenderse bien mediante dichas 
funciones y puede aplicarse a la optimización del 
diseño de sistemas ópticos con estas características. 


También es reciente la investigación de haces de luz 
con simetría cilíndrica (anulares), que generan haces 
de luz que pueden ser representados matemáticamente 
con funciones de Bessel y de ahí se han propuesto 
posibles aplicaciones [42]-[54]. A esta relativamente 
nueva rama de la óptica se le conoce como "óptica de 
haces no-difractivos” (o de haces de Bessel). De 
especial interés es la micromanipulación de partículas 
[55]-[64] donde las propiedades no difractivas de los 
llamados haces de Bessel los hacen una opción muy 
interesante para esta tecnología [65]-[88]. La 
particularidad teórica de estos haces y algunas de sus 
características han sido estudiadas en función del haz 
de luz y sus propiedades ópticas [89], debido a que 
permitirían, por así decirlo, evitar el llamado "límite de 
difracción" de algunos sistemas ópticos, lo que ha 
constituido hasta la fecha una limitante experimental 
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para varias aplicaciones tecnológicas, ya que los 
fenómenos de difracción presentes en esos casos, 
impiden mayores resoluciones. Es en este sentido, la 
importancia teórica del estudio que se presenta en 
ésta tesis adquiere relevancia, ya que actualmente es 
de gran interés, debido al potencial impacto 
tecnológico. 


En el contexto anterior, el presente trabajo de 
investigación pretende analizar con detalle los 
aspectos teóricos de la interferencia múltiple de haces 
de luz con simetría cónica, tema no tratado en la 
literatura más que muy superficialmente, y que, como 
demostraremos en los capítulos siguientes, puede 
tener importantes e interesantes aplicaciones 
prácticas. 


Los objetivos concretos de este trabajo de tesis se 
pueden resumir en los puntos siguientes: 


1. Analizar matemáticamente, y con detalle, los 
fenómenos de interferencia de varios haces 
convergentes con invariancia a rotaciones 211/N, 
mostrando las particularidades de tener esta 
simetría específica, mediante análisis en espacio 
real y espacio recíproco. 
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2. Utilizar herramientas teóricas relativamente particular la configuración propuesta en el 
novedosas en óptica, como la función de presente trabajo. 
Ambigiiedad y función de distribución de 
Wigner, que han demostrado su utilidad en otras 


ea El tratamiento matemático que se presentará en los 
ramas de la física. 


capítulos siguientes requiere de conocimientos de 
óptica física y se presupone que se conocen conceptos 
de análisis de Fourier y de Distribuciones de Espacio 
Fase. 


3. Proponer una nueva herramienta matemática, 
basada en la Transformada de Fourier y 
funciones de ambigiiedad y de Wigner, que nos 
permite obtener un espacio fase ad-hoc al 
análisis del fenómeno de interferencia, en 
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Sistemas de coordenadas cartesianas y esféricas. 
Las segundas dependen de su posición en el espacio. 
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Capítulo I 


planteamiento del problema y obtención de las 
ecuaciones de intensidad del sistema para 
polarizaciones azímutal (P,) y cenital (Es) en 
coordenadas esféricas 
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Actualmente, las investigaciones sobre 
interferencia de haces no difractivos o 
haces de Bessel e interferencia de haces 
múltiples han cobrado gran interés 
dentro de la física óptica y la 
nanotecnología, debido a su uso como 
"pinzas ópticas” en el primer caso, y al 
desarrollo de patrones de interferencia 
periódicos y cuasiperiódicos que también 
permiten la manipulación de átomos con la 
finalidad de obtener arreglos cristalinos 
o cuasicristalinos. Mediante un análisis en 
espacio real, recíproco y de fase, esta 
tesis pretende contribuir teóricamente a 
la comprensión y caracterización óptica 
de algunos de estos fenómenos. Para 
lograr esto, analizaremos teóricamente el 
caso de un conjunto de N haces de luz 
convergiendo radialmente, Figura 1.1. 
Idealmente, todos son haces de luz 
coherente con frente de onda plano 
(provenientes todos de un rayo láser) e 
incidiendo con el mismo ángulo cenital, la 
misma intensidad luminosa y fase inicial. 


Figura 1.1 Yao DL 


k, We... 
Consideramos un sistema de AY haces de luz 
convergiendo radialmente. Se muestran algunos haces 
pasando por el origen. 
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Los rayos se encuentran distribuidos uniformemente sobre el ángulo azimutal y los frentes de onda forman 


Figura 1.2 5 


Los A! haces tienen frente de onda plano, convergiendo 
radialmente en el eje Z. Esta figura está girada 180? 
alrededor del eje Z respecto de la figura anterior. Sólo se 


muestran algunos de los frentes de onda. 


pirámides de MN lados, Figura 1.2. Deseamos 
encontrar la ecuación para la intensidad del 
campo T' debida a la interferencia de los N 
haces planos en la superficie normal al eje de 
simetría. La sección transversal de los haces 
es idealmente infinita. Para encontrar dicha 
ecuación, consideremos la ecuaciones de N 
haces monocromáticos planos: 

Y lp,t)= A, expli, DARE, )] (1.1) 
donde 

p=xi+y+zk, (1.2) 

con coherencia temporal (provenientes de 
una fuente como un rayo láser), convergiendo 
radialmente, siendo 4, el brillo o intensidad 
luminosa, k, el vector de propagación y e, la 
fase inicial del enésimo rayo de luz 
respectivamente, w la frecuencia temporal y 
t la variable temporal. El subíndice n 
identifica a cada haz desde 1 hasta M. La 
función de onda total es la suma de todos los 
componentes del campo electromagnético: 
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Y.= Ny = exp|- ¡ar]Y A, expli(%, -PHE, )L, 


n=1 


y la intensidad del campo resultante es el promedio temporal del cuadrado del módulo de la amplitud compleja total: 


AS , (14) Figura 1.3 


1, bo y)=[%,88,1]- Y, ,1Jar= 


Siendo R* el espacio tridimensional. Las condiciones iniciales 
consideradas son: que el ángulo cenital es igual para todos los 
rayos y cada rayo se encuentra distribuido equidistantemente 
entre rayos adyacentes a lo largo del ángulo azimutal. La 
coherencia no será considerada en este trabajo. 
Matemáticamente, estas son las condiciones anteriormente 
mencionadas: 


O, =P Vm; (1.5) 
27m 
O, ==. 1.6 
n= (1.6) 
De la ecuación (1.4): 
NON 
m=1 n=1 


Considerando la ecuación (1.3), la ecuación (1.7) se transforma en: 


1,=Y Y A, -4,explil(£, —E,)-p+(e, —e,)l. (18) 


m=l n=1 


con lo que la variación de la intensidad resulta estacionaria. Los vectores de onda son de magnitud constante, ya que 


El vector E, —£, es paralelo al plano XY. 


estamos considerando un conjunto de haces con frente de onda plano y colineales al vector de propagación: 


17 


18 


A 277 21Ml. 13 27. |2M|. 1737, 27 
Es o sin[ó)? — en sn sin[p]; — 57 cos[ó]% , (1.9) 


donde 4 es la longitud de onda de la luz que forma el patrón de interferencia. El término k, —K_, Figura 1.3. de la 
ecuación (1.8) y con la ecuación (1.9) es igual a: 


e =k a al - 005 ) 9] Pr] sin[p]. (1.10) 


El vector k_ —k_ es perpendicular a la ondas estacionarias (perfil del patrón de interferencia) generado por el par 
de haces mm, n, paralelo al plano X-Y de registro del patrón de interferencia y exactamente a la mitad de los dos 
vectores k, y k, que lo forman. De forma que, para el producto escalar entre ambos vectores, obtenemos: 
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Polarización azimutal Polarización cenital 


Figura 1.4 A Ñ Figura 1.4 B 
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Definimos las constantes angulares para x e y 


como: Figura 1.5 A 


128 =sin|[" (rn) sip [23 lsin[ó]: (1.12) 
Úte =cos[ 0 lsin[:lsinfó]. (113) 


Consideramos dos casos de polarización: cuando el 
vector de polarización es azimutal (polarización 0 ), 
Figura 1.4A y cuando es cenital (polarización $), 
Figura 1.4B. El valor de la función de onda 
electrica es, para cada caso: 


Y (0)= Y 40, ex (a .p= == HE: ) : (1.14) 


m=1 


271ct 1 0 a ea 
Y. 1bj= A, exp ik, P2 22 +€, | .(115) Vector emésimo de propagación, eléctrico y 
magnético para polarización azimutal. 


O =-si E o, 
' =—sin sá 1 +C0S 5 Ve 
$ = 605 cos[ó); + si cos[ó]7 —sin[ól% . 


El valor de cada amplitud o vector eléctrico (4, =E,,) es Figuras 1.5 AyB: 
4,10 +=4, (=sinl6, li + cos[6, 7): (1.16) 


Siendo: 
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A, ió, j= 


5 Á, (cos[e,, ]cosló Y + sin (o. ]cos[ó]; —sin [6]%)' 
(1.17) 
El producto 4, - Á, es para el caso azimutal: 


Figura 1.5 B 


y para el cenital: 


ASA A A: í — 2sin? A o 0) 
q N 
(1.19) 


Vector emésimo de propagación, eléctrico y Sustituyendo la ecuación (1.11) y la ecuación (1.14) 
en la ecuación (1.7) obtenemos para la intensidad 


a " ” " 
magnético para polarización cenital. ) 
“9 Pp Pp del conjunto de haces con polarización azimutal: 


LIÓ= YY 4,4, l ==) in 0.0 (6.=E: )] | (1.20) 


m=1 n=1 Á 


Sustituyendo la ecuación (1.11) y la ecuación (1.15) en la ecuación (1.7) obtenemos para la intensidad con polarización 
cenital: 


1, 6)= yy AA, ( =) in os p) amd lc. -Q. y) le, —e, )] (1.21) 


m=l n=1 
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Estas son ecuaciones generales para haces con distintas amplitudes y fases. Para simplificar estas fórmulas, 
imponemos que la intensidad luminosa y la fase inicial son del mismo valor en todos los rayos (experimentalmente, 
corresponde a un conjunto de haces provenientes de un mismo haz láser, con la misma longitud de camino óptico): 
A, =A,=4:; (1.22) 
E, = €, =€. (1.23) 
La variación de la intensidad es entonces: 


LÍO )= 4 > _ 2sin ue Me. 00: »)| (124) 


m=l n=1 


LG )= 4 2 a asin Pa! cos” lle Ac. 20. ») | (1.25) 


m=l n=1 


Mediante operaciones algebraicas, se demuestra que las ecuaciones obtenidas para la intensidad son ecuaciones en el 
dominio de los reales (pues la parte imaginaria es: isin|f (0, —0 )]+isin[f(0,-6 )]=0) 


LA8|= 4? > E asi E asi (Cs =D. »)) (1.26) 


m=1 n=1 


LAG)=4? > — 2sin? ERES 1] — 2sin? EJE -0O,. »)) (1.27) 


m=1 n=1 


La ventaja de ésta formulación es que nos permite ver cada suma en términos de cuadrados de funciones, lo que es 
una formulación clásica de la intensidad. También podemos reformular como: 


1,0 V=4? eo A o (C,,x=0,. »)| (1.28) 


1,18) EE sinlóle oo] 002 los leo (C,,,x -Q.) (1.29) 


m=1 n=1 
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Las condiciones impuestas hacen al sistema invariante al girarlo un ángulo 2T1/N alrededor del eje Z para N par. 
Notemos, sin embargo, que para valores de N impar son invariantes a rotaciones r/N. Esto marca una diferencia 
entre las configuraciones de haces par e impar. Cada configuración es también invariante ante traslaciones a lo largo 
del eje Z, esto es, el plano r— 0 puede colocarse en cualquier valor de Z. Esto se debe a las condiciones ideales 
supuestas: haces de luz de frente de onda plano, fuentes en el infinito y el no considerar condiciones de frontera 
para los haces. La intensidad con polarización cenital tiene una componente horizontal y otra vertical: 


Ar 


1,467= A? cos*[9]> S cos 27m=n) cos SEO mE: x-0,.y) + 4? sin” [6] cos lc, x— O...) - (1.37) 
22 o 1 22 


Por lo que podemos identificar dos componentes fundamentales: 
1.- Paralelo al plano r—6: 
N N = 
A A O n= OQ, 5) (1.38) 
m=1 n=1 


en esta configuración, las fm componentes de las polarizaciones se encuentran distribuidas radialmente en el 
plano r—60, por lo que la llamaremos componente radial. No es propiamente la intensidad para la polarización 
azimutal, pues ésta es una configuración particular mientras que (1.38) es una componente. Podemos decir que la 
configuración azimutal sólo tiene componente radial; 


2.-  Perpendicular al plano r—6 o vertical: 


1, (Lj=40 E (a :-Q..y)| (1.39) 


m=1 n=1 


que describe un caso real sólo cuando H=77/2. Ésta la llamaremos componente axial. 


Una característica de estas configuraciones, es la escalabilidad: al variar el ángulo cenital, varía la frecuencia, 
cambiando de tamaño la configuración pero sin variar su estructura, como se muestra en la gráfica 1.1. 


23 
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Gráfica 1.1 Para las configuraciones con polarización radial y axial, al variar el ángulo cenital variamos la escala del 
patrón de interferencia, pero no su estructura. Esto es válido para cualquier N. 


y [a] EIP¿P¿1O y, 8), [42] 


618 22 50 os | SEMONSEEZMONREAA 


3 p-2.4 1.8 
3.5 p-2.8 2.1 
4 p-3.2 24 


4 


0.0804 — _, 


3.6 -32 -28 -24 -2 16 1,2 -0.8 -04 0 04 08 12 16 2 24 28 32 36 4 


A A A A A E O E E RN E E Te e EA E AMO E 


0.0641 6 -45 -4 -35 -3 25 2 15 41 05 0 05 1 15 2 25 3 35 4 45 5 
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La intensidad para la polarización /,, (5) queda, en términos de componentes radial y axial: 


1,46 )=c0s'[6]1, [+ sin*lo]7, LL), (1.40) 


y puede describirse gráficamente como: 

1, ij=c0s lol, Ali+sin*lol/, LJ. (141) 
Esta fórmula muestra las dos componentes fundamentales, 
además de hacer notar que son perpendiculares entre sí. 


Mediante las configuraciones para polarizaciones azimutal 
y cenital es posible obtener una fórmula para la intensidad 
para cualquier estado de polarización lineal entre haces, 
introduciendo un ángulo « en el plano de estas 
polarizaciones figura 1.6 y tomando dichas polarizaciones 
como componentes del vector de amplitud: 


A0, =A0 coslar]+ 44, sin[a]. (1.42) 


m m 


El producto entre dos amplitudes m, nes, entonces: 
A, Ala)= al + asi O (Go [6]sin* [a] — 0) 
Para A=MX ( M entero ) obtenemos la configuración 
: MT l 
azimutal, y para CN obtenemos la cenital. La 


intensidad para cualquier estado de polarización lineal es 
entonces: 
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Figura 1.6 3 


Mediante un vector en el plano B..0.. podemos 
generalizar las fórmulas 18) e 1%; en una para 
cualquier polarización lineal, para cada función Y. El 
ángulo «, caracteriza el estado de polarización 42, 


para dicha función 1Í2;. 
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LAB += 4 Y ( =) sn - sin” [6]sin? a) cos (AA ») | (1.43) 


m=1 n=1 


Esta configuración permite "girar" el haz de polarización alrededor del vector de propagación, mientras que la 
configuración azimutal permite “girar” el vector de propagación alrededor del vector de polarización. 


Obtendremos la configuración de interferencia para 2, 3, 4, 5 y 6 haces con P=7/2, para las configuraciones axial 
y radial. A partir de éstas haremos el análisis de la configuración azimutal y general. 
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Gráfica 1.2.2 Perfil de la intensidad para dos haces, polarizaciones radial y vertical: las intensidades están giradas 
una de otra debido a que para el primer caso, los vectores de polarización apuntan en direcciones opuestas (eje Y), y 
en el segundo en la misma dirección (eje Z). En ambos casos el máximo de intensidad es 4 A?. 


LIP¿1O0 y, 9 =35) hlP3100 y, 9 =35) 


hIP¿1 142] 
MP 311A3] 


Do yu> 
Daya 
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Gráfica 1.2.3 Perfil de la intensidad para el caso de tres haces, formando un entramado hexagonal. No son 
exactamente uno el negativo del otro como en el primer caso ya que el valor de la intensidad no es la misma; 4.5 A? 
para la configuración radial y 9 A? para la axial. 


> 2 
R[P2106 y, 9 = al L[P¿T0% y, $ = 7) 


la[P31,147] 


l[P21,002] J e 


45 


e ' 
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Gráfica 1.3.3 El patrón de interferencia sí es el mismo para ambos casos. Este es un caso de jaula óptica hexagonal. 


y, 12] IP 10, y, 0 =3),1A2] y, 12] lPJOs y, 8 =3),142] 


p  —— ——  — 


0.8 
0.4 
0.2 
—0.2 
0,4 


x, [a] 
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Gráfica 1.4.3 A Gráfica tridimensional de la intensidad para tres haces, polarización cenital. Nótese que para 


fp=1/2 se obtiene la configuración de la polarización axial. El plano horizontal que se observa corresponde al 
ángulo de extinción. 
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Gráfica 1.4.3 B Gráfica tridimensional de la intensidad para tres haces. Se observa que no hay cambio de 
frecuencia, pero sí de intensidad. Para 4=7/2 se obtiene la configuración de la polarización axial. El plano 
horizontal que se observa también corresponde al ángulo de extinción. 
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Gráfica 1.2.4 Para cuatro haces, es más evidente la diferencia entre los patrones de interferencia para las 
configuraciones radial y axial. 
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Gráfica 1.3.4 Gráfica tridimensional de la intensidad para cuatro haces, polarización cenital. Para p=7/2 se 
obtiene la configuración de la polarización axial. En este caso no hay extinción total. 
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Grafica 1.4.4 A Gráfica tridimensional de la intensidad para cuatro haces, polarización cenital. Para p=7/2 se 
obtiene la configuración de la polarización axial. En este caso no hay extinción total. 
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Grafica 1.4.4 B Cuatro haces, variando el ángulo « . 
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Grafica 1.2.5 Para cinco haces la diferencia entre los patrones se acentúa. Como se obtiene un arreglo con simetría 
10, podemos llamar al patrón un entramado de jaula óptica “cuasicristalino". 
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Grafica 1.3.5 En el origen obtenemos un anillo circular, los primeros mínimo y máximo de un haz Bessel de orden 
cero ( J, ) para la componente radial ( ángulo azimutal) y los primeros máximo y mínimo de un haz Bessel de orden 
uno ( Jo) para la componente axial. 
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Grafica 1.3.6 Variación del patrón de interferencia para cinco haces, al variar el ángulo cenital. 
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Grafica 1.3.6 Variación del patrón de interferencia para cinco haces, al variar el ángulo «. 
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Grafica 1.2.6 Para seis haces también obtenemos un arreglo hexagonal. 
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Grafica 1.3.6 Podemos ver que aunque la configuración es distinta, sin embargo los puntos centrales del entramado 
(mínimo para el caso radial y máximo para el axial) coinciden en el plano, para el mismo ángulo cenital. 
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Grafica 1.4.6 A Seis haces, polarización cenital. 


IP Xy, 0), A] IP dy, 0), A] 


E x 
: 7 
L 
de _a 
16 E 16 
a Ed , 
8 8 
az E 
16 16 
ES xr 
7 7 
q, rad Ed 4, rad Ed 
16 16 
je x 
q 3 
Ei Ei 
16 16 


CFATA UNAM 


Grafica 1.4.6 B Configuración para seis haces. 
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Grafica 1.3 La variación de la intensidad Zy 
con relación al ángulo cenital f en el plano XZ 
(espacio real), y = 0,N=(2,... ,6). La variable 
espacial x está en unidades de longitud de onda 


A 
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Serie de gráficas 1.4 Variación de la intensidad en el eje X respecto a la variación del ángulo de polarización y en 
la columna A y respecto al ángulo f en las otras dos. 
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Graficas 1.4 Plano X-Z (y=0) para N = 2, 3, 4, 5, 6. Mediante los tres vectores característicos de la función de 
onda, podemos agrupar las tres columnas según el eje de giro: 


e para la columna A, alrededor de cada vector de propagación giran sus vectores eléctrico y magnético; 
. B,alrededor de cada vector magnético giran sus vectores de propagación y eléctrico; 
e C, alrededor de cada vector eléctrico giran sus vectores magnético y de propagación. 


En los tres casos presentados, los ejes de giro respectivos están distribuidos en el plano 7-6. Cuando los ejes de 
giro se encuentran en el eje Z, todos los ejes coinciden y el patrón de interferencia sólo rota alrededor de Z. Para 
cada N, es posible obtener infinitas configuraciones intercambiando los ejes de giro. 
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La serie de gráficas de intensidad del conjunto de haces con polarización azimutal muestran líneas de interferencia 
continuas, regulares para los casos de dos, tres y cuatro haces. Para un número de haces mayor, la frecuencia se 
vuelve irregular. Las gráficas de intensidad para polarización cenital muestran extinción (ausencia de interferencia) 
y cambio de fase y bifurcaciones en la frecuencia. La extinción es total para N= 2, 3 y se debe a que los haces 
forman ángulos rectos entre ellos a un ángulo f específico: 45” para 2 haces y 54” 44' 8.197" para 3 (los tres 
vectores formando el vértice de un cubo). La extinción es parcial para N>3. El cambio de fase para la extinción 
parcial se debe a que para un ángulo /, un haz se encuentra a un ángulo recto con algunos haces, pero no con todos y 
se observa en las gráficas como bifurcaciones en la frecuencia. 


El fenómeno de extinción ocurre cuando dos haces tienen sus vectores de polarización perpendiculares entre ellos. 
Para que esto ocurra en el caso de polarización cenital, es necesario que: 


¡asi PE cos .]=0. (1.44) 


Siendo el ángulo f, el ángulo cenital al que ocurre la extinción entre dos haces m, n y es igual a: 


asi] | sj 
1 ; N 
/, =arcco = arcsin im, n=1,..., NY. (1.45) 
10, si] 20 0 si] 
N N 
Debido a la geometría del arreglo, se deben satisfacer las condiciones: 
a Md (1.46) 
2 2 
0o<6,<*. (1.47) 
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Esto se debe a que el ángulo real entre ambos vectores de polarización debe ser T/2 para que haya extinción. 
0 —60_ es mayor ya que es la proyección sobre el plano XY del ángulo real entre los dos haces considerados. La 


gráfica 1.2 muestra el número de haces contra el ángulo de extinción entre cada par de haces, para 512 haces. 
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Gráfica 1.5 A Número de 
haces N contra el ángulo de 
extinción df, entre pares de 
haces. Se grafican todas las qe _ 

combinaciones posibles. f, está E OS o 
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Para el caso de polarización azimutal, la condición de extinción es independiente del ángulo cenital $, y sólo ocurre 


l l ? Tr 37% AE 
cuando la diferencia angular azimutal es AA Notemos que esto ocurre únicamente para 


configuraciones con N haces que son múltiplos de 4. Para el caso de polarización Y, ecuación (1.43), la condición de 
extinción es: 


1-2 sn —sin” [6]sin”[ar])=0, (1.48) 


el ángulo cenital al que ocurre la extinción entre dos haces m, n es: 
. ,|Aim=n 1 
sino (m2) | 
N 2 
. -|aim—=n 
lola PR 
N 


MT , 
Para a (M entero) se obtiene el caso cenital, ecuación (1.31), y para 4 > MZ se puede expresar la 


p, 0 a la - sr quo - 1) (1.50) 


que nos define matemáticamente el caso azimutal explicado anteriormente. Podemos también definir el ángulo «,.: 


mo 


sin lin Poor! 


p,Ío,0, — 0 P=arcsin 


m 


(1.49) 


ecuación (1.49) como: 


O, =arcsin (1.51) 


51 


52 


que nos indica que los vectores de polarización de dos haces deben estar en ángulo recto para que exista extinción. 
Esta simetría entre f, y Q, es debida a que los vectores de propagación y polarización de un haz son 
perpendiculares entre sí. La gráfica 1.3 muestra el número de haces contra el ángulo de extinción entre cada par de 
haces, hasta para 512 haces, para diferentes valores de «Y. La gráfica 1.4 muestra el valor de f, respecto a la 
diferencia angular azimutal entre dos haces m, n. 
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Gráfica 1.5 B En este caso, 
a. =71/4. Aunque es la misma 
estructura de la gráfica 1.2, 
nótese que el rango de (, está 
en el intervalo de O a 1/2; esto 
significa que el ángulo de 
extinción es escalable con el 
ángulo «a. Conforme disminuye 
el valor del ángulo « al girar 
sobre el vector de propagación, 
el valor del ángulo f, aumenta, 
y para valores de a<xm/4 
éstos se salen de rango. 


UNAM 


In(P) N vs. pela =<) 


de, Pad 


0 32 64 96 120 160 192 224 256 200 920 992 5904 416 dde 400 512 


El 


9) 


54 


. 
0 
. 


A a A A E]. 
ela | a Hi ci cs o ela 
ls Aaa a. A 5 1] 
a DE = 
Es 

uN 
73 
o | 
un 

CIS E AAA q 
E 
cn 
- 

X 
! 
> 
E ¿lo 

OS AAA. E 
Z +1” 

A A A 5] 

a o o ia ¿lo 

E 
Elo 
> Elo 
E E E E E E E E E E E 

Secesjasioclesicsin so ce sjesa o ejes o 
O 
E GR onE am 8 ss 
03: > ss + D + 
ar + OY y o 3 

y Yoo + NOE 
A O 5 
E Z R gs 
Q Y < ta cc (12 0 — 
A -S9NSRagEÉE sn» 
< N m0 mn S 3 
a di MN Ñ ES 
y SN2Z 4 SES 
Sel oz dl Ye $ 
TADEO A Sa 
—302= IS a $ 

na n 
Se O 
”“"25< o Why E50NY 
sa *.o3.2 s pl 

== 
Se vuY2os 3 Rs 
$40.020 _3+ e MM 
Eure 25 3 ga 3 
Vd 32n30u yu 2D 7 


CFATA UNAM 


Cuando tenemos un ángulo de extinción para algún par de haces para la polarización cenital, al aumentar el ángulo 
azimutal también tenemos una inversión de fase para dicho par de haces. Si escribimos la intensidad total como una 
suma de intensidades parciales 


N N 
Í,= 22 Ln, , (1.52) 


la fórmula explícita del cambio de fase para cada término m, n es entonces: 


(1 — 2sin? Eo 25 C.. - 0.) - il > >m=n>N 
1, AO = A? 0, mon=2N (1.53) 


(2 sin” [pu=»)] a 1 — asin (Cs = 0.) > <m=n<N 


(1-2cos”[p]sin” 501 2sin | (cx 0..») - il 1 - 2cos' [ó]sin”[+u72]>0 
IL, ABj=A4% 0, 1 —2cos? [6]sin*[14-)]=0 . (1.54) 
(2cos [p]sin*[*4;>"]- if - 2sin” E (C,,,x— | 1 —2cos* [ó]sin*[14-2]<0 


Obtenemos una suma de senos y cosenos cuadrados, siendo los términos coseno debidos al cambio de fase. También 
tenemos términos negativos. Sin embargo, el valor total de la intensidad en cualquier punto del espacio es positivo o 
cero. La fórmula no toma valores negativos debido a los términos m= mn. Esto puede concluirse al calcular el valor de 
la intensidad en el origen, que es igual a cero, y es el mínimo minimorum de la ecuación. 
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Serie de Gráficas 1.7 Las configuraciones para 3 y 4 haces pueden utilizarse como "jaulas ópticas”, en especial las 
configuraciones radiales. Las jaulas ópticas sirven para ordenar un conjunto de átomos o moléculas en un entramado, 
en éstos casos en una cuadrícula o en un arreglo hexagonal. Los átomos se alinean en los huecos donde la intensidad 
del campo es mínima. 
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Gráficas 1.8 Cinco haces forman un anillo circular (primer mínimo y máximo de un haz Bessel), por lo con ésta 
configuración es posible obtener una "pinza óptica”. 
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Graficas 1.9 Comparando el perfil de la intensidad de un haz Bessel con perfiles de intensidad para 5 y 6 haces. El 
primer orden (máximo-mínimo) coinciden entre sí. 
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Gráficas 1.10 Para seis haces, tenemos un conjunto de haces Bessel (primer máximo y mínimo) en un arreglo 
hexagonal, formando una jaula óptica de pinzas ópticas. 
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Haciendo el número de haces muy grande dividiendo la función de onda entre el número de haces (conservando la 
misma energía), obtenemos: 


timo tim Ao cios y explik, - p]=iJ [27rsin[ó]]6 : (1.55) 
pi ES pi ARI of p|=J,[rsinló]le. (1.56) 


Esto muestra que la polarización radial genera un haz de Bessel de orden uno, y la componente de polarización axial, 
haz de Bessel de orden cero. Las intensidades correspondientes son: 


lim viÓs_ 7 ¿Dar sin[p]]: (1.57) 
mt =J [27 sinlo]]. (1.58) 

La componente cenital es una combinación de los dos haces: 
imc sn o:l2rsinlo])+ cos (917: 2rsinlo]. 159) 

y, en general, para cualquier conjunto de polarización lineal: 
a amsinloll+ Vi Resinlóll- 4 Resinlólsin [ón la] (160) 


Es posible desplazar el ángulo de extinción sobre el eje Z, girando el estado de polarización de los haces de luz. 
Variando los ángulos Y y f y el número de haces N es posible obtener configuraciones con zonas de extinción 
deseadas, para ser utilizadas como "pinzas ópticas” y "jaulas ópticas". 
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Se obtuvo una fórmula general para la intensidad de N haces de luz, con cualquier estado de polarización para dicho 
conjunto. La polarización del conjunto de haces puede por tanto, caracterizarse mediante el ángulo «. La 
información del estado de polarización se pierde parcialmente en la ecuación de la intensidad: sólo es posible 
conocer el valor Y en el primer cuadrante, no es posible determinar el cuadrante sobre el que se encuentra la 
polarización. A partir del estado de polarización, se hace una descripción del ángulo de extinción entre pares de 
haces. Mediante esta fórmula, es posible: 


Utilizando las configuraciones propuestas, obtener haces de Bessel de orden cero y uno perfectos en 
primeros órdenes, o combinaciones específicas de éstos; 


Obtener configuraciones periódicas y cuasiperiódicas variando el número de haces; 


Para una configuración con un ángulo cenital fijo, modular únicamente la amplitud de la interferencia girando el 
vector de polarización (Y sobre su eje (el vector de propagación); 


Para el caso de polarización azimutal, podemos también modular únicamente la frecuencia, girando en esta 
caso el vector de propagación sobre su eje (el vector de polarización); 


Combinando estas tres cantidades, podemos obtener configuraciones apropiadas para pinzas ópticas, jaulas ópticas, 
o como en el caso de seis haces, una jaula óptica hexagonal de pinzas ópticas. La figura 1.7 muestra un posible 
arreglo experimental para cinco haces. 
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Figura 1.7 Arreglo experimental esquemático para cinco haces. El riel permite variar el ángulo cenital, y los 
retardantes de fase variable coordinados el ángulo Y en forma sincronizada. 


Arreglo experimental para 5 haces, vista superior. 
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Capítulo IT 


ecuaciones de intensidad en espacio recíproco mediante 
Análisis de Fourier 
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El análisis de Fourier es una bien desarrollada herramienta matemática que permite descomponer una función 
periódica f,en sus frecuencias constituyentes mediante la Transformada de Fourier (TF) de dicha función 3|f]. 


Esta nos da la correspondencia de una función en espacio real x, yen su correspondiente ecuación contraparte en el 


espacio recíproco u, vy se define para dos dimensiones como: 


3..ly0oy)l= | [7(x, y)Jexpl- 2m(ux + vy)lardy. 


xX=->%0 y=-00 


Aplicando la TF a las ecuaciones (1.35) y (1.36) en función de xe y: 


S., a A |! 


m=1l n=1 


So 1,1 H]=4 YyYS, E CS, E 7 9,.9|| : 


m=1 n=1 


se obtiene el flujo luminoso: 


3. IM v=22%4 99 (1-2sin [*e-2Dójau 20, [2 +20,, ] 


m=1 n=1 


3.,[1, QU Mu, vb=28 4? yy dl -2C,,,Jó[4v +20, ,]. 


Estas dos ecuaciones también son escalables, como se muestra en la gráfica 2.1 
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(2.1) 


(2.2) 


(2.3) 


(2.4) 


(2.5) 
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Grafica 2.1 La TF de la intensidad para los casos radial y axial son escalables con el ángulo de incidencia: al variar el 
ángulo cenital, varía la escala sin variar la configuración. Al compararse con la gráfica 1.1, nótese que la escalabilidad 
varía en forma inversa entre aquella y ésta. 
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Tomando en cuenta las propiedades de la TF, con las ecuaciones (2.4), (2.5) y (1.38) obtenemos para la polarización 
cenital: 


3. 11,14]=cos"[6]S, [2,411] + sin*[o]S,,[z, 41). (2.6) 


La fórmula del flujo luminoso para cualquier estado de polarización es: 


m=1 n=1 


3.11, (PHu,v)=24 4 ze - 2(1 —sin*[6]sin*[a])sin? as -2C,, ló[2v+20,,]. (2.7) 
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Serie de gráficas 2.2 Variación del flujo luminoso en el plano recíproco U—V. Se corresponde con la serie de 
gráficas 1.4. Nota: el flujo luminoso de las dos columnas fue editado en photoshop, considerando el flujo obtenido en 
las ecuaciones (2.6) y (2.7). 


ny 
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La serie de Graficas 2.2 muestran el plano recíproco U—V para N = 2, 3, 4, 5, 6. En cada caso, los ejes de giro se 
encuentran en el plano r-6. 


e En columna A, el eje de giro es el vector de propagación; 
e B,el eje de giro es el vector magnético; 
e C,el eje de giro es el vector eléctrico. 


Nótese que para columna A, al variar el ángulo «a el flujo luminoso varía pero la frecuencia espacial permanece 
invariante (modulación de amplitud), mientras que para la columna C al variar el ángulo f varía la frecuencia pero el 
flujo luminoso permanece invariante (modulación de frecuencia). Para la columna B, el flujo y la frecuencia varían, ya 


Y Y . Y T Y . 
que ambos ángulos varían inversamente, de acuerdo con la fórmula « = pes En los tres casos, el máximo efecto se 


presenta cuando los ejes de giro respectivos se encuentran radialmente distribuidos en el plano -- 6, figura 2.1. 


La configuración con N=2"  (n=2,3,4...) polarización radial, presenta una ausencia de líneas de interferencia, 


debido a los pares de haces cuyo ángulo azimutal entre ellos es de 90”, y que, en éste caso, es independiente del 
cenital, ya que el vector de polarización siempre se encuentra en el plano gráfica 2.2-4C. La intensidad con 
polarización radial presenta un cambio de fase cuando el ángulo azimutal es mayor de 90". 


La extinción entre pares de haces ocurre cuando el ángulo entre ellos es igual a 1/2. Nótese que en todos los casos, 
ésta ocurre para u=-/24. El ángulo cenital H al que esto sucede, puede calcularse con la fórmula (1.29). 


De forma análoga a 1,16), SÍ, (61 es una combinación circular de S|1, 10! y S[Z, X2Y], ecuación (2.6). Para éste 
caso el flujo luminoso es equipotencial radialmente en el plano WU V, gráfica 2.3. En especial, la circunferencia de 


radio 2 (u? +v) es donde ocurre la extinción. 
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La topología del patrón de interferencia es distinta para los haces pares e impares. Para la configuración axial, 
gráfica 2.4 1-A los haces impares muestran A-N+1 puntos de interferencia distribuidos en la intersección de 
(N-1)/2 circunferencias con 2N radios distribuidos equidistantemente y uno más al centro del arreglo. Los haces 
pares, gráfica 2.4 1-B y C muestran A/2+1 puntos de interferencia. La configuración presenta M/2 círculos con 
2N radios, y los N puntos de interferencia de cada círculo se distribuyen alternadamente en las intersecciones, 
tanto radial como circularmente, además de uno al centro. La configuración radial presenta la misma estructura para 
ambas configuraciones, con un caso como excepción. La diferencia entre ambas configuraciones es que en la 
configuración radial los puntos se encuentran modulados en amplitud, de acuerdo al ángulo entre ellos gráfica 2.4 
2-A y B. Sin embargo, en la configuración radial con N múltiplo de 4, gráfica 2.4 2 C una circunferencia se 
encuentra exactamente en el radio de extinción, por lo que el número total de puntos de interferencia en estos 
casos son W/2-N+1. 
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Gráfica 2.3 Plano U-V para el flujo luminoso, polarización cenital. 
v.12] soe Las líneas de interferencia para cada N proyectan éste círculo de 
r intensidad. En este caso, el círculo correspondería a N= «o. La 


extinción ocurre en la circunferencia [2 (u? +v) j 


-2 -1.75 -421.25 -1-0.75-0.5-0.25 0 0.25 05 0.75 1 1252 175 2 


Serie de Gráficas 2.4 Patrón de interferencia en el plano recíproco U—V , configuraciones axial, renglón 1 y radial, 
2. La columna A muestra el patrón de interferencia para N = 9, la B, N= 10 y la columna C N= 8. Cada patrón 
muestra una configuración particular en función del número de haces y el tipo de configuración. 
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Figura 2.1 Modulación de Amplitud, frecuencia y ambas al girar dos vectores de los N campos electromagnéticos 
respectivos alrededor de un tercero. El máximo efecto se obtiene al estar los Nejes de giro distribuídos en el plano 
r—0. Si éstos se encuentran alineados en el eje Z, el giro corresponde únicamente a una rotación del sistema sin 
variaciones. 


£,, : modulación de Amplitud 


a eta 


A A O 
— A 
E,, + modulación i E,,: modulación 
de frecuencia de frecuencia y 
amplitud. 
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Capítulo III 


análisis de los sistemas estudiados en espacio fase 
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En el estudio de sistemas ópticos es muy ilustrativo representarlos simultáneamente en espacio real y espacio 
recíproco. A esta representación en ambos espacios se le conoce en óptica como espacio fase. Dentro de éstas 
representaciones, hay dos funciones útiles en especial: La función de Ambigiiedad (FA), definida en dos dimensiones 
como: 


Arno y eS | 7 +x,y + y]: 7 lx —x, y — y ]expl2om(ux” + vy")lax'dy”. (3.1) 


—o0—o0 


y la distribución de Wigner (DW), que en dos dimensiones se define como: 


Pon Wsyl= | [7er xy + 1-7 Lex, y — y lexpl- 2m(ux + vy")lax'dy”. (3.2) 


—o0—o0 


Estudiaremos el caso para dos haces incidiendo sobre una superficie, con un ángulo H, medido desde el eje Z figura 
3.1. Las formulas para los haces de luz incidiendo son: 


P.=4 exp sino — zcos|p|- 00) (3.3) 
Y, =4, AS xsinló]- zcos[p]- 0) (3.4) 
Para la polarización radial tenemos: 
A =4J: (3.5) 
A, =-Aj. (3.6) 
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Y para la polarización axial tenemos: 


A =Ak:; (3.7) 
4, =Ak. (3.8) 
Las funciones de onda son respectivamente: 
1 j ] 2 1 4 
Y, JU) =2¡4exp A (ecos[e]+er) 0 a (3.9) 
Y, [1)=24exp A (2cos[e]+er) 0 A (3.10) 
Sustituyendo (3.9) y (3.10) en (1.41) obtenemos la función de onda para la polarización cenital: 
Y, [4 =24exp 2 (cos[o] + a) ( cos|6]sin cial j+sin[6] 0 A ¿ . (3.11) 
Las intensidades para (3.9), (3.10) y (3.11) son respectivamente: 
LO =44 sin” A a +24 0 A e (3.12) 
L,(1]=44 cos” El =24* +24 00 A (3.13) 
LI4=24? -24” cos[24]cos EA . (3.14) 
La TF de la intensidad en cada caso es: 
2si 21 
3. [1,01] =246[u] 408 al aa | (3.15) 
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3,11 ()]=248[u]+ £8 
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AO e A 


3.11 14]=248[u]- 4? cos Eb al + 23 4] l 


El espacio real y el reciproco se corresponden en las variables x > u, 


Al calcular la DW y la FA para (3.9) obtenemos: 
l 2sin[g] 


"(9 09]= (a SEE +0 


Y para (3.10) se obtiene: 
4xsin[p] 


Pe.) ES (13 | =24 | 


Ao ES (13 | =14 o 


Notemos que: 


| 
ant] ta) 
| 


4 xsin|p] $ 


Wee) ES (13 | —= a ES [13] . 
La serie de gráficas 3.1 muestra las ecuaciones (3.18), (3.19), (3.20) y (3.21) para distintos valores de ángulos 


cenitales. 


y> v. 
eN 2sin[6] iS 4xsin| 6] 
| y) | y) 
24? cos En [e] +7 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 
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Serie de graficas 3.1 Gráficas de la FA y TW para distintos valores de $. El primer renglón corresponde a la 


configuración axial, y la segunda para la configuración radial. Las columnas corresponden al mismo ángulo para ambos 
casos. 
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-2 n n x.2 , =2 , a rr AN 
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Mediante la serie de gráficas 3.1, podemos generar el espacio tridimensional de la FA y WD (que para el caso axial 
es el mismo, y para el caso radial difieren por el signo). La serie de gráficas 3.2 muestra el espacio fase 
tridimensional formado con el ángulo cenital ¿ originado en (3.23). Nota: Estas gráficas fueron generadas 


parcialmente en Mathematica y editadas en Photoshop. 
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Gráfica 3.2 A El espacio 
fase generado mediante las 
TW y FA a la función de 
onda con polarización axial 
muestra la gráfica x-g4 en 
el plano u=0, y dos 

superficies 
correspondientes a la TF de 
los productos cruzados. 
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WIP. 30,U,9)], ACÉXP¿J0CU 9] 


Gráfica 3.2 B El espacio 
fase generado mediante las 
TW y FA a la función de 
onda con polarización radial 
muestra la gráfica x-g en 
el plano u=0, y dos 

superficies 
correspondientes a la TF de 
los productos cruzados. 
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Serie de gráficas 3.3 El plano x-¿ para u=0. 


6, rad a(E[P]lu = 0) 


nl ——_—_—_—_—_— —»_——_——_ 0 
1 


a 
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Tomando u=0 en (3.19) , utilizando (3.23) y comparando con (3.13) resulta: 


W 


(x,u=0) 


[P,(1)]=24 0 A =1,11)-24?. (3.23) 


Haciendo la misma operación para el caso radial en (3.18) y con (3.22)y (3.12), se obtiene: 


Wo. EXI p]=24 A 24”. (3.24) 


En ambos casos, se obtiene el término de interferencia entre los dos haces considerados. El plano x- 4 en u=0 se 
muestra en la serie de gráficas 3.3, para las configuraciones axial y radial. 


Tomando x=0 en (3.19) , utilizando (3.23) y comparando con (3.16), se obtiene: 


A í-0.1) (Y, (| SS Wecoj) ES (| E 44”8[u] +24 eE EE o as ou 200 an. S: sd Ly] (3.25) 


haciendo lo mismo para (3.18) con (3.22) y (3.15): 


a] +245(u]=3,1.(3] — (626) 


El plano u— 4 se muestra en la serie de gráficas 3.4, para x=0. 
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Serie de gráficas 3.4 Plano u—¿ para x=0. La primer gráfica corresponde al caso axial (gráfica 3.2 A) y al radial 
el segundo (gráfica 3.2 B). Obtenemos la TF de la intensidad para el caso axial. Para el caso radial, para FA y -WD 
obtenemos la TF correspondiente. 
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Utilizando las propiedades de aditividad de la WT con (3.18) y (3.20) obtenemos wm v(4): 


W(*(4| =24? a + li] + o m4] 44? cos[26] 00 A! ó[u)l, 


y con (3.19) y (3.21), A[Y(4||: 


A 


Aro [E 14 ]=44? 0 aa 24 colza] 3 + ed + dE 


En este caso, las ecuaciones no son las mismas debido a la asimetría mostrada en (3.22). 


Para el caso más general, tenemos: 


W 


As) ES [a | =44? on 


e LP. La] ]=24 a + 


47 xsin[p] ] 


2sin[6] | 
A 


+ JE - al] - 4.4? (cos? [0] + sin? [ar] cos[2o)cos 


A 


47 xsin[p] 
A 


S[u]- 24? (cos” [4] + sin” [9]cos[26)) a + e! + JE 


10 


(3.27) 


(3.28) 


(3.29) 


(3.30) 
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Capítulo IV 


hipótesis 
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La hipótesis fundamental de este trabajo se expresa matemáticamente como el uso de la siguiente función de 
distribución, que se aplicará sobre las funciones de onda estudiadas: 


Do EXC = Í Fla+x IF [x+x]exp[2rriux"] dx! = Í [FL + x]/ exp|-27.iux'] dx”. (4.1) 


Que es una TF de la intensidad, integrando sobre una variable desplazada. Aplicando o|*, | para las 
configuraciones axial y radial, obtenemos: 


3. (9, (9]=240[0]+ 4 0 A a. a + dE a]. (4.2) 
O EL] = 2408 4] +4 o au Al + a al]! (4.3) 


La serie de gráficas 4.1 muestra las ecuaciones (4.2) y (4.3) para valores específicos del ángulo cenital ¿. 
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Serie de gráficas 4.1 Planos x—u para o|*, | para distintos valores de $. El primer renglón corresponde a la 


configuración axial, y la segunda para la configuración radial. Las columnas corresponden al mismo ángulo para ambos 
casos. 


a APAP LO =01 


. xA 
5 4-3 2-1 0 1 2.3.4 8 5 4-3 2-1 0 1 2 53 


xA 
4 5 5 4-3 2-1 0 1 2.3.4 5 


7 APIO o PILI 


' 
cola 


27 AÑO =0.1 


5-4 -3-2-10 01 02.03.4050 5 44-32-1001 2 03 


4 


mn 


5 4-3 2-1 0 1 2.3.4 8 
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La serie de gráficas 4.2 A y B muestran el espacio fase tridimensional, componentes axial (4.2) y radial (4.3). Nota: 
Estas gráficas fueron generadas parcialmente en Matemática y editadas en Photoshop. 
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Gráfica 4.2 A Espacio 
Fase generado mediante la 
función 0, [Y,(13], ec. 
(4.2). La superficie 
intermedia (plana) se 
muestra parcialmente 
transparente. 
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Gráfica 4.2 B Espacio 
Fase generado mediante la 


función Me E). 
ecuación (4.3). La 
superficie intermedia 
(plana) se muestra 


parcialmente transparente. 


ña 
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La diferencia entre los planos de la serie de gráficas 4.1 y la serie 4.3 es que las primeras son planos a un ángulo 
cenital fp específico (cortes) y las segundas son la proyección x—u para todo el rango de ¿ (vista superior). Para 


28 
obtener la ecuación de ésta proyección, se utiliza la ecuación (4.2), haciendo la sustitución u = Ol 
o (Y, (1]]= 4? cos| 271 it as (4.4) 
nto) á = A" Cos | 277 x 7 u<?2, E 
con ésta ecuación se obtiene la proyección x-¿ del lado derecho de la gráfica 4.2 A. Sustituyendo u =0 en (4.2): 
Do) DES (| a 24” u= 0 (4.5) 
, aos es : 2sin[6] 
obtenemos el plano vertical de la misma gráfica. Finalmente sustituyendo u =-— ] : 
o (E, (17]= 4 cos ar | uo (4.6) 
ES a = A" cos| 27x 7 <u<0, A 


se obtiene la misma proyección x- $, pero del lado izquierdo. Sumando (4.4), (4.5) y (4.6) obtenemos el valor de la 
intensidad para la configuración axial, ecuación (3.13): 


Bo [11]. seso [9 (1), ca aa +2 00 2 O (47) 


Xx U=- x,u 
A 


Realizando la misma operación para la ecuación (4.3) se obtiene el valor de la intensidad para la configuración radial, 
ecuación (3.12): 


A] [*, ( y] DS retrata) ES 
4 A 


[aia] 


Y] =24 +24 00] 2 +] (4.8) 


x,u 


a ES ( y Es el 
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eS ; de / aos 2sin 
Para obtener la proyección x—-u para la configuración axial, se realiza la sustitución ne, en (4.4), 


de 
A =u en (4.6) y se suman ambos términos con (4.5): 


(Y, [Lp] =24 +24? cos[27xu4]  -2<u<2, (4.9) 


x,u= 


Mo [Pp + pala ) [*, (13 ]+ El sud 


correspondiendo esta ecuación a la gráfica 4.3 A. Repitiendo el procedimiento para la configuración radial: 


21 sin[$ 
A 


x,u x,u= 


CA in rato) EXA re 


ple! p=24 +24 cos[27xu+x]  -2<u<2, (4.10) 


correspondiendo esta ecuación a la gráfica 4.3 B. Para la proyección u—¿ sustituímos x=0 en (4.2) y(4.3), 
obteniendo respectivamente: 


Boo [PJ] =2408[u]4 4 as A al] =S. [1,113]: (4.11) 
Bos [PY] = 2428 [u]- 4 Al +8 A] =S, [140]. (4.12) 


En el plano x=0 obtenemos la TF de la intensidad correspondiente en ambos casos. Puesto que el rango de x es en 
principio infinito hacia ambos lados, no es posible obtener una proyección de esta perspectiva, pero sí cortes a 
distintos valores de x. Esta idea se muestra en la serie de gráficas 4.4. 
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Serie de Gráficas 4.4 Plano u—G, en el primer renglón para x=0 y en el segundo para x=54. La primer columna es 
para el caso axial, y el segundo para el radial. Para el primer renglón obtenemos la TF correspondiente, y para x=52 
los valores de la intensidad de la TF en x=512. 
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2 Ed 
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D u, pr ] D u pr ] 
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La fórmula para la configuración con polarización cenital es: 


Me Y [4 ]=24*8[u]- cos [24] 4 0 A + a + dE - al] y (4.13) 


Y para la fórmula general, obtenemos: 


D..., [9,125] =24 5[u]- 4? (cos[29]sin? [0] + cos” [a] 0 A al, A A! + o E] - (4,14) 


Aplicando la fórmula (4.1) a (1.14) y (1.15), podemos generalizar el procedimiento para Nhaces: 


Dias E (2)J-2E ase ES [6]sin? a) o (C,,,-0, 5) 8 l E rs 15 


m=1 n=1 


Que podemos sintetizar como: 


8... [%]-3É El a al aa) 


m=l n=1 
Mediante ésta función, podemos identificar las fórmulas qué ar las proyecciones o planos en distintas 
perspectivas: 


N 


Nov, r],y te] = YAA ¿cos E Ca) [E (4.17) 


m=l n=1 =1 n=1 


u—v(x=y=0) es la pane Transformada de Fourier: 


N N C,, A 20, al N N 
A A NY A, Adu a ar _ DIA (4.18) 


m=l n=1 m=1 n=1 m=1 n=1 


x—u e y—vson las superficies curvas caracterizadas por: 


DES 2) 7 YA A co a a (4.19) 


m=1l n=1 m=1 n=1 
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AOS], 20 = Ny A, A, 0 a (4.20) 


La proyección x-u de las gráficas 4.3 son invariantes a ÑN, por lo que pueden utilizarse para cualquier N > 2 y 
también para caracterizar haces de luz cónicos, entre éstos los haces de Bessel, que cuentan con la propiedad de 
ser "no difractivos”, haciendo la sustitución: 


D[P[p, os 1+ cos27r,/x* +y dul +v? + | (5.9) 
OPIP Hoc 1+ cos2r,/x* +y du? + y] (5.10) 
o[v(», )]=< 1 [te +) Jefa jv + y lu? +] (5.11) 


o[v(?, Jo 1+ (E +v ) sin*[or]- ] cos[27r.[x? + y? Ju? + v*] (5.12) 


CFATA UNAM 


La utilización de esta transformada se ejemplificará aplicándola a Y, [P,) , figura 4.1. La Función de onda es: 
4 
eS 2? - 
-24Asin[6]sin E ajo A 
—2Asin[ 6] sin in] a ap 


44 oso] cs E al y 6 A a) co E dos lo] | 


La intensidad es: 


cos 0 
; (5.13) 


Lip j= 


44* +24” cos[26|cos E sin ol +24? cos[24]cos 2 sin ol ) (5.14) 


144 co [o]eos Ho] | cos o ] 


Y su TF: 
es [2] a 
440[u]5[»]+ 
24 onl2o| [a a] +3 a qq 2] . q o ble (515) 


A E AO O 
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Aplicando (5.1) a (5.13) se obtiene: 


2 o y «2, 
O A, 
ES 


(5.13) 
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Figura 4.1 Los cuatro vectores componentes de “Y, (p,) y sus vectores eléctricos. 
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Esta configuración presenta simetría rotacional 7/4 alrededor del eje z , por lo que se obtienen las mismas gráficas 
sobre el eje x y sobre el eje y. La serie de gráficas 4.5 muestra la ecuación (5.13) para cuatro valores del ángulo 
cenital p en los planosx—u y y—v. 


La gráfica 4.6 muestra la proyección x—u, que es la misma que la proyección y—»v. 


La serie de gráficas 4.7 muestran la proyección x-4 y y-0. 


La serie de gráficas 4.8 muestra los planos u— ¿$ para x=0 y v-G4 para y=0. 
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Serie de gráficas 4.5 planos x-u e y-v para ARE Se obtiene la misma configuración para ambos 
casos. 
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Grafica u lay" ade. 
4.6 ¿ 
Proyección 
Xx—Uu e 
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Serie de gráficas 4.7 Proyecciones x-f e y-($. Este plano muestra la superposición de las interferencias 
parciales mostradas en la serie de graficas 4.5 entre los cuatro haces. La segunda gráfica muestra realzada la 
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Serie de gráficas 4.8 Planos u—¿$ y v-4 para y=0. Debido a que la superficie tridimensional de espacio fase en la 
dirección del espacio real es infinita, siempre hay que especificar en qué punto se hace el corte. La gráfica A 
muestra el plano para u—4 y v-g en x=0, y=0; la gráfica B muestra el plano u—$ en x=1, v-4 en y=1 y laC el 
plano u—f en x=5,v=-G6 en y=5. En todos los casos, la extinción ocurre para u =v= IA 


A B C 


3 up v=0 > (uv =0,x=1,y=0 > ug v=0x=5 y= 
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Los planos x-u, y-v de la serie de gráficas 4.5 
muestran la intensidad entre pares de haces en su 
frecuencia correspondiente. Por otro lado, las 
proyecciones x-f y-—0 de la serie 4.7 muestran la 
intensidad total debida a la superposición de las 
intensidades parciales. Esto es, mediante la 
transformada propuesta en este trabajo, ecuación 
(4.1) es posible obtener las intensidades parciales 
entre pares de haces, y la intensidad total resultante. 
La proyección x—u de la gráfica 4.6 es invariante a N. 
Únicamente varía el rango en el espacio recíproco u, v 
dependiendo del ángulo máximo entre pares de haces. 
Aunque no fue posible obtener el espacio 
tridimensional para ésta configuración, la serie de 
gráficas 4.9 muestra el espacio fase tridimensional 
para cuatro haces, configuraciones axial y radial. El 
espacio tridimensional correspondiente a la ecuación 
(5.13) es una combinación circular de las gráficas 4.9 A 
y B: para u, v< |V2] predomina la configuración radial, 


para u,v=|W2| se presenta una zona de extinción 


(planos verticales x-¿ó en u=+V2 e y-G en v=+4V2) 
y para u, > |V2| predomina la configuración axial. 


Finalmente, para la configuración radial, es posible 
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utilizar el mismo espacio tridimensional, sin embargo 
los planos verticales x-f e y—4G se desplazan sobre la 
variable recíproca correspondiente u ó v al girar el 
ángulo a, como se muestra en la gráfica 4.10. 


La utilización de esta nueva herramienta matemática 
nos permite comprender la interferencia como una 
superposición de interferencias individuales, siendo 
cada interferencia una superficie entre dos haces. La 
curvatura de la superficie viene dada por la TF de cada 
interferencia individual, sobre ésta superficie se 
encuentra el patrón de interferencia dependiendo del 
ángulo azimutal entre ambos haces. La superposición se 
obtiene mediante la proyección de todas las 
interferencias individuales desde la perspectiva de la 
variable en espacio real y el ángulo cenital. Desde la 
perspectiva de la variable recíproca y el ángulo cenital, 
hay que hacer cortes, ya que la superficie en el espacio 
real es infinita, y el perfil de la superficie nos da la TF 
de cada interferencia. La perspectiva desde las 
variables real y recíproca muestra la misma superficie 
para cualquier N>2, en ésta caso la superficie 
corresponde al ángulo real entre pares de haces, cuya 
configuración es la misma independientemente de su 
posición en el espacio. 


Grafica 4.9 A Espacio 
fase para cuatro haces, 
configuración radial. Los 
planos se superponen y 
generan las proyecciones 
x-Q e y-0, grafica 1.4 
4 C. Nota: Estas gráficas 
fueron generadas 
parcialmente en 
Matemática y editadas en 
Photoshop. 
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Grafica 4.9 B Espacio 
fase para cuatro haces, 


configuración radial. 
Nota: Estas gráficas 
fueron generadas 
parcialmente en 


Matemática y editadas en 
Photoshop. 


Gráfica 4.10 (página 116) 
Proyección de la variación 
del plano de extinción al 
girar el ángulo « . 
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Capítulo V 


Discusión, perspectivas y conclusión 
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El trabajo de tesis comprende la caracterización 
interferométrica de N haces de luz convergentes bajo 
ciertas condiciones: coherentes entre sí, con 
invariancia a traslaciones en el eje Z y con la misma 
amplitud y fase. Inicialmente se consideraron dos 
casos: funciones de onda con polarización paralela al 


plano r—8, azimutal (P,) y con polarización cenital 


(»,) y partir de éstas dos, se dedujo la existencia de 


una componente axial de polarización (1) y otra radial 


) y (1) son fundamentales 


ra 


). Las configuraciones ( 


en el sentido de que no son composición de ninguna 
otra: por el contrario, la función de onda con 
polarización E es una composición “circular” de las 


(4). 


configuración (P,) solo tiene componente (0). 


mientras 


configuraciones (1) y que la 


Las 


componentes axiales tienen sus componentes de 
polarización alineada sobre el eje Z en la misma 
dirección hacia arriba. Por otro lado, la componente 
radial presenta los vectores de polarización 
distribuídos apuntando radialmente hacia afuera. Es 
debido a estas configuraciones que la componente axial 
presenta un máximo en el origen, mientras que la radial 
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un mínimo. Para dos y tres haces, los patrones de 
interferencia presentan la misma estructura, pero 
están invertidos: Donde (1) presenta máximos, (1) 


tiene mínimos y viceversa. Siendo (P,) una composición 


de las dos fundamentales 
)cos* [6] +P(1)sin*[6], para dos y tres 


circular 
en)" 


haces hay un punto en que ambos patrones de 
interferencia se anulan, presentándose una extinción 
total. Como el máximo de intensidad es el mismo para 
ambas configuraciones con dos haces (447), la 
extinción en éste caso se presenta para $, = 45". Como 
la intensidad máxima para tres haces es (4.54?) para 
la configuración radial y (94?) para la axial, la 


extinción total se presenta para /, = 54" 44' 8.197". En 


los demás casos, las configuraciones axial y radial no 
presentan la misma estructura, por lo que sólo hay 
extinciones parciales. 


Utilizando las configuraciones (P,) y (P,) como 


componentes, se obtuvo una función de onda (P,) más 
general, mediante la introducción un ángulo «a entre 


las configuraciones (P,) y (P,) que permite "girar" el 


vector de polarización. Cada configuración permite la 
modulación del patrón de interferencia. La 
configuración (P,) tiene sus vectores eléctricos 


distribuidos radialmente, sobre cada uno de éstos 
giran los vectores de propagación y magnéticos para 
modular la frecuencia. 


Uno de los objetivos del presente trabajo fue el de 
analizar las características únicas que la simetría 
cilíndrica (haces convergentes) pudiese tener, 
obteniéndose las ecuaciones que describen el patrón 
de interferencia en el plano de r-6, así como sus 
transformadas de Fourier, que las describen en el 


espacio recíproco u, v. Para el caso LAB describimos 


las zonas de extinción entre pares de haces, cuando 
estos se encuentran en ángulos rectos entre ellos en 
términos de sus componentes angulares axial y radial. 
Esto ocurre para todos los casos de Na diferente (, 
el ángulo cenital al que ocurre la extinción entre dos 
haces lo llamamos 6,. Para los mismos pares de haces, 
también encontramos que el patrón de interferencia es 
diferente para ángulos menores a 90” que para ángulos 
mayores: específicamente, hay un cambio de fase igual 
a tt. Esto es evidente en configuraciones de dos y tres 
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haces, en donde hay un ángulo que presenta extinción 
total. Para la configuración con polarización (P,) esto 


ocurre cuando el par de haces se encuentra a 90% 
sobre el plano, ocurriendo en configuraciones donde el 
húmero de haces Nes múltiplo de 4. 


Por medio de la aplicación de la transformada de 
Fourier, notamos que la extinción total en (E se debe 


a que en estas configuraciones hay una sola línea de 
interferencia en las gráficas. Para N > 3 hay más de 
dos, y ninguna coincide con otra en la interferencia 
nula al mismo ángulo ff, debido a que la zona de 
extinción ocurre a un valor invariante para todos los 


casos, Vu? +v? =v/24” ya que depende únicamente de 
la longitud de onda. Sin embargo, el cambio de fase 
ocurre para toda M. El cambio de fase se representa 
como un húmero negativo. 


Se aplicaron las Transformadas de Distribución de 
Wigner y Función de Ambigiiedad a la función de onda 
de N haces, obteniéndose una representación 
tridimensional en espacio fase. Para dos haces, esta 
representación muestra la superposición de ondas 
mediante el plano x-f en u=0, y dos superficies 
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curvas que muestran el valor de le frecuencia para 
cualquier ángulo cenital. 


También se obtuvo otra representación matemática 
para la intensidad D[f(x, y)K(x, y), (u, v)?, basada en 


las funciones de Fourier, distribución de Ambigiiedad y 
de Wigner, que son funciones intermedias (como se 
llama también a las funciones en espacio fase), y que 
constituye la hipótesis fundamental del presente 
trabajo. La función propuesta puede equipararse a una 
Transformada de Fourier, aplicada sobre una variable 
de integración desplazada. Es importante recalcar que 
la función de distribución se aplica a la ecuación de 
onda de los haces de luz que la componen. De manera 
similar a la representación del espacio fase 
tridimensional para FA y WD, esta función proporciona 
información sobre el espacio real y el recíproco 
simultáneamente en función del ángulo cenital. La 
diferencia radica en que las gráficas nos permiten 
representar a la función de distribución como un 
conjunto de superficies, en la que cada superficie 
representa un ángulo azimutal O entre dos haces. La 
variación de la frecuencia entre pare de haces se 
puede observar en cada superficie en función del 
ángulo cenital f. Mediante su utilización para N=4, 
se ejemplifica cómo el espacio fase tridimensional 
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las 


x—Uu 


presenta información sobre intensidades 
individuales mediante planos para 
distintos valores de (; las proyecciones x-4 e y-G 
muestran la composición total de la intensidad, y los 
planos u— $ en x=0 y v-ó en y=0 muestran el perfil 
de la TF de la intensidad. Las proyecciones x-u e 
y—v muestran la invariancia entre el espacio real y el 
recíproco para distintos valores de M_ pues esta 
superficie (o su parte correspondiente) es la misma 
para cualquier valor de N mayor de 2, y representa el 
ángulo real entre dos funciones de onda cuyos vectores 
de propagación coinciden en un punto. 


e y-v 


Estas fórmulas, en general, permiten obtener gráficas 
utilizando herramientas comunes de computación, en 
específico el programa Mathemática, y los paquetes de 
graficación que contiene. En todos los casos, se utilizó 
una alta resolución para las gráficas. 


El trabajo teórico realizado nos ha permitido 
caracterizar el comportamiento interferométrico de N 
haces de luz con invariancia a giros 211/M. Este estudio 
abre nuevas posibilidades de investigación y 
desarrollos tecnológicos: 


Diseño de “jaulas ópticas" (optical cages) 
periódicas con haces de luz con N = 3, 4, 6 y 
cuasiperiódicas con N= 5, y > 7; 


Caracterización de haces cónicos de luz 
convergentes: N > 5 nos da un anillo de luz para el 
primer máximo/mínimo central, y es posible utilizarlos 
como "pinzas ópticas” (optical tweezers). 
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El concepto de modulación de la frecuencia, amplitud o 
ambas mediante el giro de cada uno de los tres 
vectores característicos de la función de onda en el 
plano r—-8 permite manipular más eficientemente el 
patrón de interferencia. 
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